Spontaneous vortex phase in ErNi 2 11 B 2 C 
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Below Twfm ~ 2.75K, a weak ferromagnetism and superconductivity coexist in ErNi2B2C. In 
the present paper, a possibility of a spontaneous vortex phase was examined directly by small 
angle neutron scattering technique. An external magnetic field was applied to align ferromagnetic 
domains so that the internal magnetic field induced by the weak ferromagnetism becomes uniform 
and uniaxial over the bulk sample. When the field is removed, the flux line lattice subsisted below 
Twfm but disappeared above it. A pinning effect is also discussed. We conclude that the spontaneous 
vortex state is realized in ErNi2B2C. 
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Ferromagnetism (FM) and superconductivity (SC) are 
both cooperative quantum phenomena which occur at 
q = 0, but a widely accepted consensus is that they 
are mutually exclusive. The pioneering theory on this 
issue by Ginzburg Q stimulated a number of theoreti- 
cal works which follow continuously since 1950's up to 
now. At present, the coexistence of such two states is 
believed to occur only when an internal magnetic field 
TJint induced by a spontaneous FM moment is smaller 
than the critical magnetic field H c for a type I SC and 
H C 2 for a type II SC, respectively ||. Moreover, for a 
type II SC, an appearance of a spontaneous vortex phase 
(SVP) is predicted if H int further satisfies a condition, 
H c i < Hi n t < H C 2 By contrast, intensive exper- 

imental studies on FM SCs were activated in the late 
1970's by the discovery of two families of magnetic in- 
termetallic SCs, RE-Rh 4 B 4 and RE-Mo 6 X 8 (RE=rare 
earth) Some compounds in these two families were, 
indeed, confirmed to show both FM and SC. When the 
FM ordering was established in these systems, however, 
SC was immediately destroyed, and the coexistence of 
these two states was observed only in a very limited nar- 
row range of temperature like AT ~ 0.1K just below Tq 



Recently, the study of magnetic SCs has been re- 
vived due to the discovery of two new magnetic SCs. 
One is an intermetallic SC, RETM 2 B 2 C (TM=Ni, Pd) 
UlUn]] , while the other is a family of Rutheno-cuprate 
compounds, RuSr 2 RE 2 Cu 2 Oio and RuSr 2 RECu 2 8 
^2| , ^3|Jl4| . By magnetization, specific heat and/or resis- 
tivity measurements, the coexistence of weak ferromag- 
netism (WFM) and SC was suggested for ErNi 2 B 2 C |j| 
as well as Rutheno-cuprates |Il],|l2],|l3],[L4| , and a possibil- 
ity of the SVP on these compounds was theoretically dis- 
cussed [M. A couple of experimental reports concluded 



the existence of the SVP in Rutheno-cuprates p7| , |i"8| ], 
but such conclusion should be taken with caution be- 
cause it was based on only macroscopic measurements. In 
fact, very recent neutron scattering studies proved that 
RuSr 2 RECu 2 0s (RE = Gd, Y) has actually an antifer- 
romagnetic ground state, and under an applied field, the 
field-induced WFM order coexists with SC p9 20|. From 
this lesson, we would like to emphasize that it is cru- 
cial to confirm coexistence of FM and SC by means of 
microscopic measurements. 

ErNi 2 B 2 C is the first material which is microscopically 
confirmed to exhibit the coexistence of WFM and SC be- 
low Twfm ~ 2.75 K [fHH). Thc results of detailed anal- 
yses of the magnetic structure in the WFM phase will 
be reported elsewhere [p3[ . Note that the previous neu- 
tron scattering studies revealed that localized magnetic 
moments on Er atoms first form a transversely-polarized 
planar-sinusoidal (SDW) structure below Tn ~ 6 K with 
an incommensurate propagating vector q\ ~ 0.553a* 
and an amplitude of the modulated moment /i = 7.8/^b 
p^ , p5[ . The crystal structure of ErNi 2 B 2 C belongs to 
tetragonal I4/mmm symmetry at room temperature, but 
it transforms to an orthorhombic phase below Tn p6[ . 
The magnetic moments in the ordered state are parallel 
to the orthorhombic 6-axis [ p7[ . 

To further examine the influence of the WFM ordering 
to SC, we have studied temperature (T) dependence of 
lower critical field H c \ in details [p2[ . As explained in 
the inset of Fig. 1, three different magnetization curves, 
( (1) vir : virgin curves, (2) dec : H-decreasing curves, 
and (3) inc : H-increasing curves) were measured with 
applying an external magnetic field H cxt either parallel 
to the b (|| 6) or c (|| c) axes, and the corresponding H CI 
were determined for each curves as the field where the 
two guided lines are crossing as exemplified in the in- 
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FIG. 1. (color)Temperature dependence of lower critical 
field: (a)H^ t , and (b)H^ xt , respectively. The inset of the 
top panel shows a typical example of M — H curves observed 
at T = 4 K for After Ref. [22]. 



set. The difference between H c \ and H CI is very small. 
Therefore, one can approximetly regard H CT as H c i- (For 
detailed explanations, see Ref. (2gj.) It should be noted 
here that a magnetization curve of the magnetic SC is in- 
evitably affected by H lnt produced by the magnetic com- 
ponents. For H cxt < H c i, however, the influence of iTnt 
is expected to be very weak because the spins which con- 
tribute H lnt are limited to those which reside within the 
penetration depth A from the surface, thereby justifies 
the following analyses and discussions. 

In Fig. 1 is depicted thus obtained T dependence of 
H cr . The upper and lower panels correspond to behavior 



of H CI for H^} xt and H^ xt , respectively [22( ]. All six curves 
exhibit a weak dip near Tn- They are attributed to the 
magnetic SDW and a structural phase transition at Tn- 
By contrast, the drastic influence was observed at Twfm. 
H CI determined by M v - lr curves with H^. shows a precip- 
itous decrease, whereas all other curves show an abrupt 
increase. By a field-cycling procedure, we confirmed that 
these behaviors are completely reproducible. Accord- 
ingly, the drastic changes of H cr below Twfm must be 
intrinsic to ErNi2B2C. 

The mechanism of the decrease of H cr below Tvvfm can 
be interpreted as follows. Remember that the decrease 
of H cr was observed only in the Mdec curve for H^ xt , but 

He 

not observed for H^ xt . This distinct difference can be 
understood by taking accounts of the domain effect and 
of the magnetic anisotropy of ErNi2B2C simultaneously. 
ErNi2B2C has a strong basal plane magnetic anisotropy 



and, as aforementioned, it undergoes the tetragonal-to- 
orthorhombic structural transition at Tn- Hence, after 
the zero field cooled (ZFC) procedure, there are four-fold 
magnetic and crystallographic domains below Tn within 
the basal plane. Below Twfm, the WFM order naturally 
forms a domain structure, and the internal field H int pro- 
duced by the spontaneous FM component would be can- 
celed out to zero. It explains why the system shows a 
complete Meissner effect in all M v ; r curves. In another 
words, after the ZFC procedure and for H cxt < H c i, the 
system behaves just like it does not know the existence of 
the WFM ordering. When H cxt is further increased and 
exceed the H c ±, it starts to penetrate into the sample. 
Here, if the H ext is parallel to the 6-axis, it starts to align 
not only magnetic domains but also crystallographic do- 
mains. Once FM domains are aligned, H mt remains finite 
and subsists with keeping parallel to H ext until H ext de- 
creases to T. The finite H lnt supplements H cxt , and H c \ 
(H CI in the present study) is drastically decreased below 
Twfm- Because of a strong in-plane anisotropy, however, 
.ffixt cannot align magnetic domains, and a decrease of 
-ff cr is not observed. 

We interpreted that the increase of H cr below Twfm 
indicates a possibility that the WFM enhances the SC 
Meissner phase. Since the WFM coexists with SC, it 
may not be intuitively unreasonable to suppose that the 
formation of the WFM microscopically enhances a sta- 
bility of the SC phase. If the WFM favors the SC, a 
new interaction must be invoked between localized FM 
moments and conduction electrons pc| . 

The most important result here is that, at lower T's, 
H CI determined by the Afd GC curves with H^ xt seems 
to go to practically zero with decreasing T. Since the 
vortex phase is expected to exist in the field range for 
H cxt > H c i, the present result strongly suggests a possi- 
bility that, when the H ext is reduced from high field, the 
vortex state remains even H cxt is reduced to T, which 
should be regarded as SVP by definition. In what fol- 
lows, we will report the direct observation of the flux line 
lattice (FLL) by small angle neutron scattering (SANS) 
technique. 

We used the same Boron- 11 substituted single crystals 
of ErNi2 11 B2C which were used in our previous studies 
pl| , [22"| ]. Note that the T c , Tn and Twfm of our sample 
were determined to be 8.6 K, 6.0 K and ~ 2.75 K, re- 
spectively |29| , |30[ ] . The remnant magnetization at 1.8 K 
was ~ 0.34^s/Er. SANS experiments were performed 
by using the SANS-U spectrometer at the Cl-2 port in 
the JRR-3M guide hall, JAERI, Tokai, Japan. The B 
and Er have high neutron absorption factor. In order to 
reduce the absorption of neutrons, we substitute natu- 
ral B with U B isotope(99 at.%), but we could not use 
Er isotope. Therefore, to enhance the transmission of 
cold neutrons, we needed to slice the large single crystals 
into pieces with dimensions of 4 x 1 x 10mm 3 . Incident 



2 




FIG. 2. (color) FLL diffraction patterns at 1.6 K (left col- 
umn) and 4.0 K (right column) observed for ZFCOT, 0.5T, 
0.1T and FDOT procedures from the top to bottom panels, 
respectively. Here, ZFC and FD denote " zero-field-cooled" 
and "field-decreased", respectively. A background data taken 
at T = 10 K was subtracted. The crystalline axes are shown 
in the left-top panel. 

neutrons with wavelength A n = 11 A were selected by a 
velocity selector and detector position after sample was 
16 m. The sample was mounted in the cryostat with 
a horizontal magnet with a field direction and incident 
beam to be parallel to the sample b axis. The data are 
collected by a 128 x 128 pixelated 2D position sensitive 
detector. Because of low intensity, it took 12 hours to 
obtain one FLL pattern. First the sample was cooled in 
zero field (ZFC). In order to align the magnetic domains, 
H cxt of 2.8 T is applied parallel to the b axis and it is 
removed. The FLL measurements are repeated at 1.6 K 
below Twfm ~ 2.75 K and at 4.0 K above Twfm, respec- 
tively. Our idea is, if SVP is realized in the WFM phase, 
then we should be able to observe the FLL signal at 1.6K 
but not at 4.0 K when H cxt is reduced to T. 



FIG. 3. Azimuthal dependence of FLL intensity. Data be- 
tween 22 < r < 90 pixel are used for the integration, (a) 
T = 1.6 K and (b) T = 4.0 K, respectively. 

Figure 2 shows typical FLL diffraction patterns in 
ErNi^E^C. As a background the 10 K data in the para- 
magnetic normal state was subtracted To omit the 
strong background scattering, the data in a small angle 
region (in radius range of r < 20 pixel) are masked in 
the figures. The sample crystallographic orientation is 
indicated in the left-top panel. At T after the ZFC 
procedure (ZFCOT), shown in the top 2 panels, only FM 
diffuse scattering is observed along the a* axis in a small 
q region. After applying H cxt as large as 2.8 T, it was re- 
duced down to 0.5 T, and then to 0.1 T (middle 4 panels). 
The diffuse scattering observed at ZFCOT is strongly sup- 
pressed at 1.6 K but not at 4.0 K, demonstrating that 
magnetic domains are aligned by ff G xt at 1.6 K. In addi- 
tion, the FLL signal is clearly observed in all four cases. 
We noticed that the q position for the FLL signal at 1.6 
K is larger than that at 4.0 K, indicating that the system 
feels an extra magnetic field at 1.6 K. This is consistent 
with our idea that £Zi nt by WFM supplements H cxt and 
affects the position of vortices. The FLL patterns after 
field decreasing to T (FDOT) are depicted in the bot- 
tom 2 panels. At 1.6 K the FLL signal is still visible but 
it is missing at 4.0 K. 

To further demonstrate the FLL signals, a portion of 
azimuthal dependence of scattering intensity is depicted 
in Fig. 3. Data in a radius range between 22 < r < 90 
pixel are used for this integration and 9 is an angle from 
the a* axis. For ZFCOT, only the diffuse components 
along the a* axis were observed at around 8 ~ deg. 
At 0.5 T and 0.1 T, however, the diffuse component was 
drastically reduced at 1.6 K, while subsists at 4.0 K. In 



addition to the diffuse components, clear peaks are ob- 
served at around 9 ~ ±40 deg, indicating the existence 
of FLL signals. Furthermore, for FDOT, the FLL peaks 
remain around 9 ~ ±40 deg at 1.6 K but completely dis- 
appear at 4.0 K, indicating that the vortices are swept 
out from the sample at 4.0 K at least until the FLL sig- 
nals disappear from the experimental q range. 

Here we need to discuss the pining effect. Recently, 
Dewhurst et al. pointed out that the bulk pining effect 
becomes important below Tn by using a miniature lo- 
cal Hall array |Q. Furthermore, Gammel et al. have 
demonstrated that a pinning force increases rapidly in 
the WFM phase by transport and magnetization mea- 
surements p3|. These results indicate a possibility that 
the enhanced pining mechanism below Twfm may trap 
vortices, but the comparatively weak pining force above 
Twfm would release them, hence the FLL signals were 
observed only at 1.6 K (below Twfm) but not at 4.0 
K (above Twfm)- They argued, however, that the en- 
hancement of the pining force was observed not only in 
the in-plane direction but also along the c axis, and they 
were of comparable order. This behavior has an impor- 
tant implication to interpret our results. Namely, in our 
study, the decrease of H CT was observed only with H^ t 
but not with H^, exhibiting a clear anisotropy. This 
result strongly indicates that the decrease of H CI is due 
to the existence of if; n t but not due to the pining effect. 

To give further confidence, one may expect that a similar 
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experiment with H£ Kt would be valuable. The previous 
FLL measurements with this field configuration, however, 
showed that the FLL rotate away from the applied mag- 
netic field direction below Twfm 1^1 , which implys that 

He 

it is impossible to observe the FLL with a small H^. 
Based on these considerations, we conclude that the pin- 
ning effect is not essential to the present case, and that 
the SVP is realized in ErNi 2 B 2 C. 

Finally in Fig. 4 is depicted schematical arrangements 
of vortices in the real space at 1.6 K for the FDOT pro- 
cedure. The field which corresponds to the observed q 
position is calculated to be about 1600 G. The FLL is 
oriented with respect to the nearest neighbor direction 
parallel to the crystalline a axis, being consistent with 
the previous work at 4 K with a lower field by Eskildsen 
et al. p5[ . From the values of H c i (H CI in the present 
study) and iJ C 2 for our sample, the coherence length £ 
and the penetration depth A at 2 K are roughly estimated 
to be 200 A and 900 A, respectively. Here to calculate A, 
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for H l £t and the M, 



curve, because for 
, the observed H cr is already re- 



we employed H c 
the case of H^ t and Md 
duced by H- mt from the field the system feels. Reflecting 
lower T c of our sample, the calculated values are slightly 
larger than those reported in references |36|,|37],|38|] . Fur- 
thermore, our detailed analyses on the magnetic struc- 
ture in the WFM phase revealed that the FM moment 
appears on the (100) plane with the periodicity of about 




FIG. 4. Schematical picture of the vortex arrangement in 
the real space at 1.6 K in the SVP state for the FDOT proce- 
dure. 

35 A p3| . The calculated A is about 26 times larger than 
the periodicity of the FM sheets. The relation between 
these length scales is also depicted in the figure. 

In summary, the vortex state in the WFM phase of 
ErNi2B2C was directly examined by the SANS technique. 
With a single FM domain, we observed FLL even at 
H cx t—0, and conclude that the SVP is indeed realized 
in ErNi 2 B 2 C. 
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